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ABSTRACT 

The Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) comprises deep multi- 
colour (u* g'r'i' z') photometry spanning 154 square degrees, with accurate photometric red- 
shifts and shape measurements. We demonstrate that the redshift probability distribution func- 
tion summed over galaxies provides an accurate representation of the galaxy redshift distribu- 
tion accounting for random and catastrophic errors for galaxies with best fitting photometric 
redshifts z p < 1.3. 

We present cosmological constraints using tomographic weak gravitational lensing by 
large-scale structure. We use two broad redshift bins 0.5 < z p ^ 0.85 and 0.85 < z p ^ 1.3 
free of intrinsic alignment contamination, and measure the shear correlation function on an- 
gular scales in the range ^1—40 arcmin. We show that the problematic redshift scaling of 
the shear signal, found in previous CFHTLS data analyses, does not afflict the CFHTLenS 
data. For a fiat ACDM model and a fixed matter density O m = 0.27, we find the normal- 
isation of the matter power spectrum erg = 0.771 ± 0.041. When combined with cosmic 
microwave background data (WMAP7), baryon acoustic oscillation data (BOSS), and a prior 
on the Hubble constant from the HST distance ladder, we find that CFHTLenS improves the 
precision of the fully marginalised parameter estimates by an average factor of 1.5 — 2. Com- 
bining our results with the above cosmological probes, we find ft m = 0.2762 ± 0.0074 and 
cr 8 = 0.802 ± 0.013. 

Key words: galaxies: distances and redshifts - galaxies: photometry - techniques: photomet- 
ric - methods: analytical - large-scale structure of Universe 
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1 INTRODUCTION 

Weak gravitational lensing by large-scale structure provides valu- 
able cosmological information that can be obtained by analysing 
the apparent shapes of distant g alaxies that have been co h erentl y 
distorted by foreground mass (Bartelman n & Schneider! 1200 lb . 
Since weak lensing is sensitive to the distance-redshift relation and 
the time-dependent growth of structure, it is a p articularly use- 
ful tool for constraining models of dark energy dAlbrecht et all 
|2006|; iPeacock et ail 120061 : lAlbrecht etail 120091) . To measure the 
contribution of dark energy over time, the lensing signal must be 
measured at several redshifts , thi s is known as weak le nsing to- 
mography (see for example, iHul Il999t iHutereii |2002[) . Several 
observations of weak lensing tomography have been complete d 
jBacon et all 120051 : ISemboloni etail |2006t iMassev et all 120071). 
Most r ecently a study of the COSMOS survey bvlSchrabbac k et alj 
d2010h found evidence for the accelerated expansion of the Uni- 
verse from weak lensing tomography. 

Redshift information is vital to weak lensing interpretation 
since the distortion of light bundles is a geometric effect and 
the growth of structure is redshift-dependent. Weak lensing data 
sets necessitate the use of photometric redshifts due to the large 
number of galaxies they contain. Spectroscopic redshifts typically 
exist for a small and relatively-bright fraction of galaxies, pro- 
viding a training set for photometric redshifts at brighter mag- 
nitudes. Several approaches for determining the redshift distri- 
bution of galaxies have been used in past weak lensing stud- 
ies. Many early s tudies (see for example. IVan Waerbeke et al. , 
20021: Ijarvis et all 120031: lHamana et all I2OO3I; iBacon et all 12003; 
Van Waerbeke et all 120051; Ifioekstra et all 120061 : iBeniamin et all 
20071 ; iFu et all [2008), lacking multi-band photometry, relied on 



external photometric redshift samples such as the Hubble Deep 
Field North and South, and the CFHTLS-Deep fields. Due to the 
small area of these fields, sampling variance was an important, but 
often neglected, source of error in these studies, as presented by 
IVan Waerbeke et all J2OO6I) . 

Current and planned weak lensing surveys have multi-band 
photometry enabling photometric redshift estimates for all galax- 
ies. Methods for measuring photometric redshifts use various 
model-fitting techniques with the goal of finding a match between 
the observed photometry and template galaxy spectra which are 
displaced in redshift and convolved with the optical response of 
the filter set, telescope, and camera. Depending on the set of 
photometric filters, degeneracies can exist between different tem- 
plate spectra at different redshifts. We refer to large errors in the 
best-fitting parameters due to mismatches under these degenera- 
cies as catastrophic errors. The effect of catastrophic errors on 
weak lensing paramet er constraints h a s been investigated in severa l 
stud ies, for example iMa et alj d2006l) ; iBernstein & Hutererl d201Cfr ) 
andlHearin et alj fold) . Using a detailed Fisher matrix analysis, 
iHearin etal 1 20 id) show the importance of properly characterising 
catastrophic errors to dark energy parameter constraints u sing weak 
lensing tomography. The implication of neg lecting these errors is 
not well known, although Ipfearin et alj d2010h argue that there are 
many factors governing the final impact on dark energy parame- 
ters and each survey needs to be carefully considered to make any 
definitive statement. It is clear that catastrophic errors will become 
increasingly important in the next generation of weak lensing cos- 
mic shear surveys. 



ysis of the Canada-France-Hawaii Tel escope Lensing Survey! 
(CFHTLenS), with redshifts measured in lHildebrandt et aL 1 2012h 
using the Bayesian photometric redshift code (BPZ. iBemtezl 
2000). The BPZ analysis of the CFHTLenS photometry uses a set 
of 6 recalibrated s pectral energy distribution galaxy templates from 
Capak et al. (2004) and a magnitude de pendent prior on the redshift 
distribution (see Hilde brandt et alll2012l for further details). If the 
galaxy template set and prior used are an accurate and complete 
representation of the true galaxy population at all redshifts, then 
the probability distribution function (PDF) calculated using BPZ 
determines the true error distribution. The redshift distribution of 
a galaxy sample can then be calculated from the sum of the PDFs 
to determine an accurate redshift distribution that includes the ef- 
fects of both statistical and catastrophic errors. This is in contrast 
to the standard method of using a histogram of photometric red- 
shifts taken from the maximum of the posterior. We test the accu- 
racy of the summed PDFs with overlapping spectroscopic redshifts 
at bright magnitu des and with resampled COSMOS-30 redshifts 
dllbert et all l2009) at faint magnitudes. In both cases we also assess 
the level of contamination between redshift bins u sing an angular 
cross-correlation technique (Benjamin et al .,2010). 

Previous CFHT Legacy Survey (CFHTLS) results were found 
to be biased, underestimating the shear at high redshifts, thus re- 
quiring the addition o f a nuisance parameter when model fitting 
(Kilbi ngeretalll2009h . Furthermore, field selection discriminated 
based on cosmology dependent criterion possibly resulting in con- 
firmation bias. The CFHTLenS catalogues we use in this paper 
have been thoroughly tested for systematic errors. These tests are 
cosmology insensitive and were completed without any cosmo- 
logical analysis of the data (Hevma ns et all l2012bh . One of the 
primary goals of this paper is to demonstrate that the redshift 
scaling of the shear is consistent with expectations. We limit our 
cosmic shear analysis to two broad redshift bins in order to ob- 
tain parameter const raints that do not depend on the model l ing of 
intrinsic alignment dHeavens et ail |2000| : ICroft & Metzlerl l200d; 
Iffirata & S eliak, 2004]). A study of cosmological constraints from 
CFHTLenS with several redshift bins, account i ng for intrinsic 
alignm ent is presented in iHevmans et alj d2012ah . iKilbinger et alj 
J2012I) present a thorough investigation of 2D cosmic shear, in- 
cluding a comparison of all popular second order shear statistics. 
ISimpson et all d2012h use the tomographic shear signal presented 
in this paper to constrain deviations from General Relativity on cos- 
mological scales. 

CFHTLenS has an effective area of 154 square degrees with 
deep photometry in five broad bands u* g'r'i'z' and a 5a point 
source limiting magnitude in the z'-band of i' AB ^25.5. These data 
were obtained as pa rt of the CFHTLS, whic h completed observa- 
tions in early 2009. IHevmans et al. I d2012bl) present an overview 
of the CFHTLenS analysis pipeli ne summarizing th e weak lens- 
ing data processing with THELI dErben et alll2012h , shear mea- 
surement with lensfit (Miller et all l2012h and ph otometric redshift 
meas urement from PSF-matched photometry dHildebrandt et all 
I2OI2I) using BPZ. Each galaxy in the CFHTLenS catalogue has 
a shear measurement e b s , an inverse variance weight w, a PDF 
giving the poserior probability as a function of redshift, and a 
photometric redshift estimate from the pe ak of the PDF z D . The 
shear calibration corrections described in lMilleretal.l(l2012l) and 
IHevmans et alj d2012bl) are applied and we limit our analysis to the 
129 of 171 pointings that have been verified as having no signif- 



In this paper we present a tomographic weak lensing anal- 
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icant systematic errors thro ugh a series of cosmo logy-insensitive 
systematic tests described in lHevmans et alj d2012br) . 

This paper is organized as follows, in Section [2] we use a se- 
ries of tests to determine whether the PDFs are sufficiently accu- 
rate to determine the redshift distributions for the many different 
science analyses of the CFHTLenS data set, and then apply our 
findings to the first tomographic analysis of the CFHTLenS data 
set in Section [3] We investigate the effect of non-linear modelling 
of the mass power spectrum and baryons on our tomographic weak 
lensing results in Section [4] Section |5]contains our concluding re- 
marks. 



2 GALAXY REDSHIFT DISTRIBUTIONS DETERMINED 
FROM THE PHOTOMETRIC REDSHIFT PDF 

When considering the redshift of an individual galaxy, a best-fitting 
redshift must be measured from the PDF, typically corresponding 
to the peak of the PDF. If many galaxies are considered, the sum 
of their PDFs can be used as an estimate of the redshift distribution 
instead of the distribution of best-fitting redshifts. We show in this 
section that, by using information from the entire PDF, we achieve 
an accurate model of the redshift distribution. The accuracy of the 
PDFs is not known a priori since this depends on whether the tem- 
plate spectral energy distributions and prior information are a rep- 
resentative and complete description of the galaxies in the survey. 

We compare the summed PDFs against several other meth- 
ods of measuring the redshift distribution. These methods include a 
comparison with the overlapping VVDS and DEEP2 spectroscopic 
redshifts (see Section l2~ll . statistical resampling of the CFHTLenS 
photometric redshifts using the COSMOS-30 redshifts (see Sec- 
tion [Z2](> an d a photometric redshift contamination analysis (see 
SectionEfl. 

We divide the data into six redshift bins and measure the red- 
shift distribution of each. We are limited in the total number of 
bins by the pairwise contamination analysis, which breaks down for 
larger numbers of bins (see Section |231 for a more detailed discus- 
sion). The redshifts are most reliable in the range 0.1 < z p < 1.3 
where comparison to spectroscopic redshifts, for i' < 24.5, shows 
the scatter to be 0.03 < <7Az < 0.06, w ith an outlier rate of less 
than 10 per cent ( Hildebrandt ~etaill20 1 2h . Here o\ z is the variance 
in the value of Az, which is given by 



Az = 



(1) 



where z p and z a are the photometric and spectroscopic redshifts, 
respectively. 

The redshift bins are chosen such that each bin is approx- 
imately four times wider than the photometric redshift error 
0.04(1 + z p ). This is done to avoid excessive contamination be- 
tween adjacent bins. For z p > 1.3 there are only a small number of 
galaxies, making a subdivision of this range difficult. The six bins 
are as follows: 



Bin 1: 


0.00 < 


Zl 


s; 0.17 


Bin 2: 


0.17 < 


Zl 


< 0.38 


Bin 3: 


0.38 < 


Z3 


< 0.62 


Bin 4: 


0.62 < 


Z4 


sj 0.90 


Bin 5: 


0.90 < 


25 


< 1.30 


Bin 6: 


1.30 < 


«6 





2.1 Comparison with spectroscopic redshifts 

We begin by investigating the redshift distribution given by spectro- 
scopic redshifts. Spec troscopic redsh i fts fro m the the VIMOS VLT 
Deep Survey ( VVDS.lLe Fevre et allfcoOSl) and the DEEP2 galaxy 
redshift survey dNewman et aUl2012l) overlap with CFHTLenS and 
were used to test the photometric redshifts. For a given photomet- 
ric redshift bin we can select those galaxies that have spectroscopic 
redshifts and examine their redshift distribution. The spectroscopic 
sample is complete for i < 22.0, dropping to ~90 per cent com- 
pleteness for i' < 23.0. We adopt the latter cut to ensure that 
there are a sufficient number of galaxies for our analysis. The cata- 
logues are also cut to exclude objects on masked regions and those 
that are flagged as stars. Stars are selected with star_f lag (see 
lErben et all l20ll for more details). Due to the dithering pattern, 
which ensures that exposures exist between individual CCD chips, 
there is a variable number of exposures over a single pointing. This 
changing photometric depth is difficult to account for when con- 
structing a random catalogue with the same properties, which is 
necessary for the contamination analysis presented in Section l2~3l 
To avoid this complexity, a final cut is made to select galaxies on 
areas of the sky that were detected during every exposure, and ran- 
dom objects are placed only in these areas. We do not expect this 
to bias our results as there is no correlation between the physical 
properties of a galaxy and which part of the CCD mosaic it was 
observed on. 

A comparison of redshift distributions for i' < 23.0 is pre- 
sented in Figure[T] For each redshift bin we show the redshift distri- 
bution predicted by the summed PDF (solid line), and the spectro- 
scopic redshift distribution (dashed line). The PDFs of all galaxies 
within a given redshift bin are summed and the resulting distribu- 
tion normalised to obtain the solid line. If the summed PDF is a 
good representation of the true error distribution, then we would 
expect this distribution to agree with the redshift distribution mea- 
sured with the spectroscopic redshifts. 

We use the Kolmogorov-Smirnov two-sample test (KS test) to 
determine if the two distributions in Figure [T] are consistent with 
being drawn from the same population (det ails of this test can be 
found in, for example, IWall & Jenkinsl l2003h . Before performing 
the test we adopt as a discriminating criterion a significance level 
of a = 0.05. The P-values found from the KS test are presented 
as labels in Figure [T] We find that the distributions for the first five 
redshift bins are consistent with having been drawn from the same 
population at a significance level of a = 0.05. However, we can re- 
ject the null hypothesis for the last redshift bin at the same level of 
significance, indicating that the two distributions are significantly 
different. This is indicative of the large uncertainties in the photo- 
metric redshifts at z D > 1.3 and confirms the choice of this cut-off 



made bv lHildebrandt et all 12012). 



2.2 Comparison with COSMOS photometric redshifts 

The agreement at bright magnitudes shown in Figure[T]is encourag- 
ing; however, the majority of lensing studies include fainter galax- 
ies, for example a magnitude limit of i' < 24.7 is adopted fo r 
the measurement of CFHTLenS galaxy shapes dMilleretal.Ll2012l) . 
Therefore, we wish to investigate the redshift distribution with this 
deeper magnitude cut. The spectroscopic redshift sample described 
in Section |2~T1 cannot be used for comparison since the complete- 
ness of this sample drops sharply beyond i! ~ 23. Instead, we 
use th e COSMOS-30 photometric redshift catalogue ( Ilbe rt et all 
2009), which is accurate to (ta z — 0.012 due to 30 bands of wave- 
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Figure 1. Comparison of the predicted redshift distributions within each broad redshift bin, labelled z;. A magnitude cut of i' < 23.0 is used for comparison 
with spectroscopic redshifts. Solid lines (pink) show the summed PDFs for all galaxies within a given redshift bin. Dashed lines (green) show the spectroscopic 
redshift distribution. The listed P- values are the result of a two-sample Kolmogorov-Smirnov test of the distributions, we adopt a significance level of a = 0.05 
rejecting the null hypothesis that the two distributions are drawn from the same population for the highest redshift bin. 




Figure 2. Comparison of the predicted redshift distributions with a magnitude cut of i' < 24.7. Solid lines (pink) show the summed PDFs for all galaxies 
within a given redshift bin. Dot-dashed histogram (cyan) shows the result of resampling the CFHTLenS redshifts using the constructed conditional probability 
P(z 3 o\z p ). The P-values are the result of a KS test, we reject the null hypothesis for the highest redshift bin at a = 0.05. 



length coverage from the ultraviolet to the mid-infrared. Quoting 
values for the Subaru i-band, the COSMOS-30 data are 99.8 per- 
cent complete for i < 2 5.5, and have a 5a point source limiting 
magnitude of i ~ 26.2 dllbert et all 12009). A resampling proce- 
dure is used to estimate the redshift distribution of deep CFHTLenS 
galaxies based on the distribution of COSMOS-30 redshifts with 
CFHTLS overlap. 

Although the 1.6 square degree COSMOS field contains the 
one square degree CFHTLS-Deep field D2, there are no overlap- 



ping CFHTLS-Wide fields. Therefore, it is not possible to directly 
match CFHTLenS galaxies to objects in the COSMOS-30 cata- 
logue. This issue can be circumvented in a novel way using th e 
photometric catalogue of D2 provided in Hildeb randt et al .1 (2009). 
Using the fact that the photometric systems for the CFHTLS- 
Wide and Deep data are identical we add random Gaussian noise, 
scaled to simulate data taken at CFHTLS-Wide depth, to the multi- 
colour magnitude estimates in the D2 photometric catalogue. Using 
artificially-degraded catalogues generated in this way, we calculate 
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a Wide-like photometric redshift estimate z v using the maximum o f 
the posterior distribution as described in (IHildebrandt et alj,l2012h . 
This is done for each D2 object in a catalo gue matched to the 
COSMOS-30 catalogue of lllbert et al] J2009t) . employing an as- 
sociation radius of 1.0 arcsec. The COSMOS-30 redshifts in the 
matched catalogue are labelled z 30 . 

This matched catalogue of z 30 and noise-degraded Wide-like 
z p estimates can then be used to acquire information about the joint 
probability distribution of COSMOS-30 and CFHTLenS redshifts. 
We generate 100 realisations of the artificially degraded Wide-like 
ca talogues, running the Ba yesian photometric redshift estimation 
of Hildebrandt etld] d2012l) for each realisation. Using this ensem- 
ble of (230, z p ) pairs, we construct a two-dimensional histogram of 
galaxy number counts in square bins of width 0.0025 in redshift 
for both z 30 and z p . This histogram is then used as an empirical 
estimate of the conditional probability density function P(z 30 \z p ) 
and allows us to estimate the corresponding cumulative probabil- 
ity distribution function P(< z 30 \z p ) for each z p bin. Then, us- 
ing inversion sampling from P(< z 3 q\z p ) with a uniform pseudo- 
random number generator, samples of redshifts distributed accord- 
ing to P(z 30 \z p ) can be drawn. 

With the assumption that P(zt Tue \z p ) = P(z 30 \z p ), the con- 
tamination in tomographic redshift bins can be estimated by resam- 
pling CFHTLenS redshifts according to P(z 30 \z p ). The resulting 
redshift distributions predicted from this method are given as dot- 
dash lines in Figuref2] and the summed PDFs are presented as solid 
lines. Note that the fine structure seen in the resampled redshifts is 
due to structures in the COSMOS field and does not represent real 
structures in the distribution of CFHTLenS galaxies. The small size 
of COSMOS means that it is limited by sample variance and indi- 
vidual clusters are able to leave an imprint in the resampled galax- 
ies. This only affects the fine details of the resampled redshifts lead- 
ing to a breakdown of the assumption P(ztiue\z p ) = -P(z 3 o|z P ) for 
small redshift intervals. 

We again adopt the null hypothesis that the two distributions 
are drawn from the same population. Using a KS test we find that 
the null hypothesis can be rejected at a significance level of a = 
0.05 for the z p > 1.3 redshift bin, but not for any of the other bins. 
The P- values found from the KS test are presented as labels in Fig- 
ure [2] Our results again confirm that the CFHTLenS photometric 
redshifts of the z p > 1.3 galaxies are unreliable. However, we find 
no evidence that the galaxies at z p < 0. 1, which is the lower lim it 
for the high-confidence redshift range (Hild ebrandt et all l2012h , 
are unreliable. This is likely because our lowest redshift bin ex- 
tends to z p = 0.17 and is therefore dominated by galaxies with 
well-measured photometric redshifts. 



2.3 Redshift contamination from angular correlation 
functions 

In order to further test the accuracy of the photometric redshift 
PDFs, we measure the reds hift contaminat i on us ing an angular 
cross-correlation technique dBeniamin et al 1 l2010h . This method 
has few assumptions and is sensitive to any contamination between 
redshift bins. Since it only relies on the angular correlation func- 
tion of the galaxies it is independent of the other methods used and 
serves as a critical test. 



2.3.1 Overview of method 

Galaxies cluster in over-dense regions, leading to an excess in the 
number of pairs found at a separation 8 when compared to a random 



distribution of points. The two-point angular correlation function 
oj(6) quant ifies this excess probab ility of finding pairs. A common 



quantifies tnis excess probabilit; 
estimator (Landv & Szalavl[T993h is 

(AA)fl N * N R _ (Afl)a N K _ (DjR)e Nr 
(RR)e NiNj (RR)e N (RR)e Nj 



(2) 



where (D[Dj)e is the number of pairs separated by a distance 8 
between data sets i and j, (RR)e is the number of pairs separated 
by a distance 8 for a random set of points, (D[R)e is the number 
of pairs separated by a distance 8 between data set i and a random 
set of points, TVr is the number of points in the random sample, 
and Ni (Nj) is the number of points in data sample i (j). The auto- 
correlation is described by the case i=j, and the cross-correlation 
by the case i^j. Our analysis would hold for any estimator of the 
angular correlation function. 

In the absence of magnification, galaxies in non-overlapping 
redshift bins should not be clustered with one another. Therefore, 
clustering between these bins should be consistent with a random 
distribution of points, resulting in uij =0. Adjacent redshift bins 
will have a small positive wy owing to galaxy clustering at their 
shared edge, which becomes more pronounced for narrow redshift 
bins. If any non-zero angular cross-correlation is detected between 
the photometric redshift bins, they must share galaxies with similar 
redshifts. 

As shown in Benjamin et al.l J20 !('■). this simple realisation can 
be exploited to estimate contamination between photometric red- 
shift bins. The reader is referred to that work for the full details of 
the method. Here we present only a few key equations and concepts 
before we apply the method to the CFHTLenS data. 

The contamination fraction, /y is defined as the number of 
galaxies contaminating bin j from bin i as a fraction of the total 
number of galaxies N?~ which have a spectroscopic redshift that 
lies within redshift bin i. If there is no overlap or contamination 
between redshift bins, /ij = when i 7^ j, and N? — N°, where 
N° is the total number of galaxies which have a photometric red- 
shift that lies within redshift bin i. In the standard case of overlap- 
ping photometric redshift bins, the contamination fraction relates 
the observed number of galaxies in each photometric redshift bin 
N° to the true underlying number of galaxies JVj as follows; 



fN?\ 

N% 



(hi 

fl2 



hi 

f-22 



/ml 
/m2 



\/lm /2m 



,/mm J 



(3) 



where m is the number of redshift bins and fa = 1 — / ik - 
We determine the contamination fractions /ij from measurements 
of 0J°j, the observed two-point correlation function between photo- 
metric redshift bins i an j. For two redshift bins it can be shown 
that, 



(jj) M| - M+<jk (jf) /2i(i - A2) 

(l-/l2)(l-/2l) + /l 2 /21 ' 



(4) 



When considering more than two redshift bins we measure the con- 
tamination fractions for each pair of bins in turn. This pairwise ap- 
proximation assumes that higher-order contamination can be safely 
ignored, that is, the angular cross-correlation is not affected by the 
mutual contamination of the pair of bins by another redshift bin. 
As the number of bins increases or if the contamination fractions 
become large, this assumption is no longer valid and the method 
breaks down. 

Once the contamination fractions /y have been measured, we 
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can invert the contamination matrix in Equation [3] to determine 
the true underlying number of galaxies in each redshift range 
from the observed number of galaxies in each photometric redshift 
bin N°. The true redshift distribution n l (zj) for each photometric 
redshift bin i is then calculated over the full redshift range, sampled 
at each redshift Zj from 

n i (z j )=f ii N?. (5) 
2.3.2 Contamination analysis 

For each pointing the angular correlation function is measured us- 
ing the publicly available code ATHEN^E ATHENA employs a tree 
data structure to increase the speed of pair counting at the cost 
of accuracy. The level of approximation is parametrized by the 
opening angle. Larger values indicate larger approximation^ with 
an opening angle of zero representing no approximation. Tests of 
ATHENA against a more simplistic and robust algorithm are used to 
determine that with an opening angle of 0.03 we are making at most 
a one per cent error on the angular correlation function. This value 
of opening angle is used when measuring the angular correlation 
function. 

We measure the angular correlation function for in six angular 
bins spaced logarithmically on the range 0.15 < 8 < 30 arcmin. 
Above 30 arcmin the signal is very small providing little additional 
information. For each pointing the contamination fractions are es- 
timated via the angular correlation function as outlined in Sec- 
tion [234] The covariance is estimated via a bootstrap technique, 
with an additional contribution coming from the field-to-field vari- 
ance for the angular cross-correlations. The details of the maxi- 
mum likelihoo d technique and covar iance matrix are presented in 
Appendix A of lBeniamin et alj ( l2010h . The likelihoods for the con- 
tamination fractions from each field are then combined with equal 
weighting. 

The following matrices contain the measured contamination 
fractions with 68 per cent confidence regions. All values are mul- 
tiplied by one hundred for ease of viewing. For the bright sample, 
i' < 23.0, we find, 



/ii = 



/65±5 


4± 1 


< 1 


< 1 


< 1 


6±6 




28 ±4 


87 ±3 


8 ± 2 


< 1 


< 1 


7± 7 


1 ± 1 


7±2 


85 ± 3 


9±2 


1±1 


7±7 




1 ± 1 


< 1 


6 ± 2 


85 ±2 


38 ± 6 


5±5 




< 1 


< 1 


< 1 


4± 1 


56 ±7 


29 ± 12 




\3±1 


< 1 


< 1 


< 1 


3±1 


18 ± 16/ 


a cut of i! 


< 24.7 


we measure, 








/42±8 


3±2 


< 1 


1±1 


< 1 


9±6\ 




32 ±4 


75 ±3 


8 ± 2 


< 1 


< 1 


1 ± 1 




< 1 


18 ±2 


79 ± 3 


7± 2 


< 1 


1 ± 1 




4±4 


< 1 


11 ± 2 


78 ±3 


20 ±3 


3±3 




< 1 


< 1 


< 1 


9± 2 


73 ±4 


42 ±5 




\18±4 


2±1 


< 1 


4± 1 


5±3 


36 ±6/ 





(6) 



(7) 

Many of the contamination fractions are one per cent deviations 

2 Since we expect the non-diagonal contamination fractions to be small the 
matrix should be diagonally dominant and therefore invertible. 

3 http://www2.iap.fr/users/kilbinge/athena/ 

4 Galaxies are grouped together into nodes in the tree data structure based 
on angular position. The structure is a hierarchy with the nodes on top con- 
taining more galaxies. The opening angle determines when to descend to 
lower nodes and higher spacial resolution. 



from zero, which is expected since we have this level of uncer- 
tainty in our estimation of the angular correlation functions. Note 
that the i th column contains the location of all bin i galaxies. Due 
to the pairwise treatment of redshift bins columns do not sum to ex- 
actly 100 per cent. These matrices are extremely well conditioned 
with condition numbers of 8.72 and 6.35 for the i! < 23.0 and 
il < 24.7 cases respectively, indicating that matrix inversion is nu- 
merically stable and does not contribute a significant uncertainty to 
the solution of equation[3) 

With the contamination fractions measured, the true number 
of galaxies in each redshift bin can be calculated from equation [3] 
The redshift distribution is then found from equation|5] This is done 
with a Monte Carlo procedure for finding global solutio ns to the 
contamination matrix presented in Benjami n et alj J2010l) . We can 
now compare our contamination results with those found from the 
spectroscopic redshifts and the COSMOS-30 photometric redshifts. 
However, since our contamination results exist in only six redshift 
bins we must also sum the distributions shown in Figures Q] and [2] 
within these six redshift bins. 

For i' < 23.0 we present our contamination results in Fig- 
ure [3] For each redshift bin we show the redshift distribution pre- 
dicted by the PDFs (crosses with solid lines), the contamination 
analysis (filled circles with dotted lines), and the spectroscopic red- 
shift distribution (filled squares with dashed lines). The horizontal 
error bars on all points denote the width of the redshift bins . The 
vertical scale is proportional to the number of galaxies but uses 
arbitrary units. For the summed PDF, the vertical error bar is calcu- 
lated as the standard deviation of the summed PDFs for 1000 boot- 
straps of the galaxies within each bin Z[. Note that given the large 
number of galaxies in each bin (~40, 000 — 600, 000) the statisti- 
cal error of the summed PDF is very small. The vertical error bars 
on the contamination results enclose the 68 per cent confidence re- 
gion which comes from a procedure for fi nding global solutions to 
the contamination (Benjamin et al. 1 120101) . The vertical error on the 
spectroscopic redshift distribution is taken as the standard deviation 
from 1000 bootstraps of the spectroscopic catalogue. For both cases 
where bootstraps are used we verified that 1000 bootstraps yields 
stable error estimates. For each bootstrap, objects are sampled with 
replacement and the resulting redshift distributions measured, the 
total number of galaxies sampled is equal to the number in the orig- 
inal catalogues. 

Figure [3] shows the predicted redshift distribution for each of 
the six redshift bins used. The contamination points for a given 
sub-plot are contained within the corresponding row of the con- 
tamination matrix in equation [6] For example, the top row shows 
that the the majority of galaxies from bin 1 remain in bin 1, with 
fu = 65 ± 5 per cent. Contamination from other bins is less than 
the per cent level except for the neighbouring bin /21 = 4 ± 1 per 
cent and the highest redshift bin /si = 6 ± 6. Keep in mind that 
the relative heights of points in the z\ sub-plot do not follow these 
contamination values because the contamination /ij represents the 
number of galaxies in bin j from bin i divided by the true number 
in bin i. However, investigating the matrix in relation to Figure [3] 
can help in grasping the presented information. We expect that the 
spectroscopic redshift distribution is the true distribution, assum- 
ing that the limited area of the spectroscopic samples does not bias 
the results, which is a reasonable assumption for our purposes. The 
contamination model is in poor agreement with the spectroscopic 
sample in the z\ and zs sub-plots. For z\ the contamination model 
underpredicts the contamination from bin 2 to bin 1, underestimat- 
ing /21 as evidenced by the discrepancy between the contamination 
point and the spectroscopic point in the second bin of the z\ sub- 
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Figure 3. Comparison of the predicted true redshift distribution within each broad redshift bin, labelled zj. A magnitude cut of i' < 23.0 is used for 
comparison with spectroscopic redshifts. All horizontal error bars denote the width of the redshift bin and points are offset horizontally for clarity. Crosses 
with solid lines (pink) denote the summed PDFs when integrated within a given broad redshift bin, the error is calculated as the standard deviation from 1000 
bootstrap samples. Filled circles with dotted lines (blue) show the result from our contamination analysis with 68 per cent confidence region. Filled squares 
with dashed lines (green) show the spectroscopic redshift data integrated within each broad redshift bin. The error is the standard deviation of 1000 bootstrap 
samples. 



: 1 1 1 1 1 1 1 1 1 1 1 

; z 1= (o.oo,o.i7] 


z 2 =(0.1 7,0.38] 


i 1 1 i 1 1 1 1 i i i 

z 3 =(0.38,0.62] 


resampled z 
i • ! contamination 


f-I— »i 




i — * — i summed PDFs 


i-m 





♦ 






■hi 1 






a-. { :m . 


EPS 








1 i i | i i ii | | | 

z 4 =(0.62,0.90] 


1 i i | i i i i | | | ■ 

z 5 =(0.90,1.30] 


1 i i | i i i i | | | 

z 6 =(1 .30,7.00] 




P— ■ 






1— X— 1 




• 
















• 




1 ■ < 

f 




i — ■ I 




m. 1-9 1 






I ■ ■ I ■ ... I I I 




_ 



0.1 



0.4 



0.9 
Z 



1.6 



2.5 3.5 



0.1 



0.4 



0.9 
Z 



1.6 



2.5 3.5 



0.1 



0.4 



0.9 
Z 



1.6 



2.5 3.5 



Figure 4. Same as Figure fJ]except for the following differences. A magnitude cut of i' < 24.7 is used. Filled squares with dot-dashed line (cyan) show the 
resampled COSMOS-30 data integrated within each broad redshift bin. The error is given as the standard deviation of the 100 low-resolution reconstructions 
(see Section l2~2V 



plot. Similarly the contamination is oveipredicted for /12 which 
is seen in the first bin of the Z2 sub-plot. This represents a fun- 
damental degeneracy in the angular cross-correlation method. Al- 
though an angular cross-correlation is detected between these two 
bins, unless the angular auto-correlations have significantly dif- 
ferent slopes, the method cannot distinguish easily between bin 1 
galaxies contaminating bin 2 or vice versa. A similar degeneracy 
explains the discrepancies in the z& sub-plot, there we see that fse 



predicted by the angular cross-correlation is too low and fes in the 
Z5 sub-plot is too high. The contamination between these bins is 
detected but the direction of scatter is misidentified. 

We use a KS test to determine if the distributions in Figure [3] 
are consistent with being drawn from the same population. Since 
there are three distributions and the KS test is a two-sample test we 
apply it to each pair of distributions. Furthermore, due to the small 
number of bins we must rely on tabulated critical values which exist 
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for very few significance levels, therefore we are not able to list 
the P-values for each redshift bin. For each pair of distributions 
we find that we cannot reject the null hypothesis (drawn from the 
same population) at a significance level of a — 0.05 for any of the 
redshift bins. 

We present the results of the contamination analysis for i' < 
24.7 in Figure [4] The summed PDF and contamination results 
are presented similarly to Figure [3] The resampled redshift s using 
COSMOS-30 are given as the dot-dash line and filled squares. The 
vertical error on the COSMOS-30 points is taken as the standard 
deviation of the 100 low-resolution resamplings (see Section [2~2l l. 
If we compare the contamination results to the resampled redshifts 
the greatest discrepancies are for the zi, zs, and ze, sub-plots. The 
fi2 — /21 and /s6 — /bs degeneracies noted for the bright sample 
above appear again in Figure |4] Additionally the resampled red- 
shifts predict a larger /ei and smaller fie than do the contamina- 
tion results which can be seen in the first and last bins of the Z\ 
and zs sub-plots respectively. The contamination analysis predicts 
that a significantly lower number of galaxies belong in bin 5 com- 
pared to the other methods, see bin 5 of the 25 sub-plot. However, 
this is not due to scattering of bin 5 galaxies elsewhere, note that 
/ss = 73 ± 4 per cent, instead the contamination analysis simply 
predicts fewer galaxies occupying this bin. To determine if these 
differences are statistically significant we use a KS test. We again 
have three distributions and apply the test between each pair. For 
each pair of distributions we find that we cannot reject the null hy- 
pothesis (drawn from the same population) at a significance level 
of a = 0.05 for any of the redshift bins. 

When using the finely binned spectroscopic and resampled 
redshifts in Sections 12.11 and 12.21 we were able to reject the null 
hypothesis for the high redshift bin z p > 1.3. The smallest P-value 
found was in the high redshift bin when comparing the summed 
PDF with the spectroscopic distribution in Figure[T] Performing the 
same comparison with these distributions when summed within the 
six redshift bins of the contamination analysis we are not able reject 
the null hypothesis. The coarse binning required by the contamina- 
tion analysis has reduced the statistical power of the test. However, 
the contamination analysis provides a complementary estimation of 
the redshift distribution, which agrees well with the other estimates 
and strengthens our confidence in the summed PDF as an accurate 
measure of the redshift distribution. 

We conclude that the summed PDF can be used to estimate 
the redshift distribution for t he high confidence red shift range 
0.1 < z p < 1.3 determined by Hilde brandt et al1d2012h . The com- 
parison with the resampled COSMO-30 redshifts and the contam- 
ination analysis for i' < 24.7 suggest that an accurate estimate 
of the redshift distribution, including statistical and catastrophic er- 
rors, can be obtained from the sum of the PDFs. T his result suggests 
that th e model galaxy spectra and priors used in lHildebrandt et aU 
(2012) are a fair and sufficiently complete representation for the 
population of galaxies studied here. 



3 WEAK LENSING TOMOGRAPHY 

In this paper we present an analysis of the CFHTLenS tomographic 
weak lensing signal using two broad redshift bins and compare our 
results with a 2D analysis over the same redshift range. Setting our 
analysis in a flat ACDM cosmology framework, the initial aim is to 
use the consistent results we find between successive tomographic 
bins as a demonstration that the CFHTLenS catalogues are not sub- 
ject to the redshift-dependent systematic biases that were uncov- 



ered in an earlier analysis of CFHTLS data (Kilbi ngeretaT1,l2009h . 
This cosmology-dependent demonstration is the last in an extensive 
series of tests, which investigate the robustness and accuracy of the 
CFHTLenS catalogues. We stress, however, that this analysis was 
performed after t he conclusion of a series of cosmology-insensitive 
tests presented in lHevmans et alj d2012bl) and the photometric red- 
shift accuracy analysis presented in Section [2] Most importantly, 
no feedback loop existed between this cosmology-dependent test 
and the systematics and image simulation tests that determined the 
calibration corrections and the subset of reliable data that we use 
from the survey. 

We choose to use two broad mid-to-high redshift bins for our 
tomographic analysis in order to reduce the potential contamina- 
tion to the signal from intrinsic galaxy alignmen ts (see, for ex- 
ample, He avens et aUl2000l ; lHevmans et alj . l2012al . and references 
therein). We estimate the expected contamination of the measured 
weak lens ing signal using t he line ar tidal field intri nsic alignment 
model of Iffirata & Seliakl d2004l) . and following iBridle & Kingl 
( 2007) b y fixing its amplitud e to the observational constraints ob- 
tained bv lBrown et alj d2002h . By limiting the redshift bins to pho- 
tometric redshifts 0.5 < z p ^ 0.85 and 0.85 < z p ^ 1.3, we esti- 
mate that any contamination from intrinsic alignments is expected 
to be no more than a few per cent for each redshift bin combina- 
tion. We therefore ignore any contributions from intrinsic align- 
ments in this analysis as they are expected to be small in compari- 
son to our statistical errors. Note that a low level of contamination 
would not be expected if we instead used the 6 narrow redshift bins 
that were analysed in the redshift contamination analysis in Sec- 
tion [2jiVe_r£ejient_a_fine 6-bin tomographic analysis of the data in 
Hevmans e t~al] d2012ah where the impact of intrinsic galaxy align- 
ments is mitigated via the simultaneous fit of a cos mological model 
and an intrinsic alignment model. The findings of iHevmans et alj 
(2012a) support the approach taken in this paper to neglect the con- 
tribution of intrinsic alignments for our choice of redshift bins. 

The 2D lensing analysis presented here is restricted to the 
same redshift range used in our tomographic analysis 0.5 < z p ^ 
1.3. We measure the shear correlation function on angular scales 
from ~1 to ~40 arcmin. The upper limit is set by our ability to 
measure the covariance matrix from simulations, see Section [3.3.1| 



3.1 Overview of tomographic weak lensing theory 

The complex weak lensing shear 7 = 71 + i"/2 , which is directly 
analogous to the complex galaxy ellipticity, can be decomposed 
into two components: the tangential shear 7 t and the cross com- 
ponent 7 X . These are defined relative to the separation vector for 
each pair of galaxies, with jt describing elongation and compres- 
sion of the ellipticity along the separation vector and 7 X describing 
elongation and compression along a direction rotated 45° from the 
separation vector. The following shear-shear correlation functions 
can then be computed: 

(8) 

where galaxy pairs labelled i,j are separated by angular distance 
■& — — i9j |. If 1? falls in the angular bin given by 8 then Ay=l, 
otherwise Aij=0. The labels k,l identify redshift bins. The summa- 
tion is performed for all galaxies i in bin k and all galaxies j in bin 1. 
The contribution of each galaxy pair is weighted by its inverse vari- 
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ance weight WiWj. This gives greater significance to galaxy pairs 

with well-measured shapes. 

Shear calibration is perform ed as described in iMiller et"aH 

d2012h and lHevmans et al] ( l2012bh . This signal-to-noise (S/N) and 
size-dependent calibration includes an additive (c) and a multiplica- 
tive (m) correction term as follows; 

7° bs = (1 + m)j 



+ c. 



(9) 



An average additive correction of 2 x 10 is found for 72. The 
additive correction for 71 is found to be consistent with zero. The 
multiplicative correction to £± is found by calcula ting the weighted 
correlation function of 1 + m (M iller et aUl2012l) . 



1 + K (tt) = . 



(10) 



The shear correlation functions £± are then corrected by dividing 
them by 1 + K. 

The shear-shear correlations can also be expressed as filtered 
functions of the convergence power spectra 



d££J 0/4 (£6)P^\£), 



(11) 



where J n is the n 1 order Bessel function of the first kind and 
£ is the modulus of the two-dimensional wave vector. These can 
be related to line-of-sight integrals of the three-dimensional matter 
power spectrum 



p k,i m _ 9 h n n 

W- 4c 4 



a (x) \JK\X) 



(12) 



where c is the speed of light, fi m is the matter energy density, Hq 
is the Hubble constant, /if (x) is me comoving angular diameter 
distance out to a distance x, Xh I s me comov i n g horizon distance, 
o(x) is the scale factor, and Pg is the 3-dimensional mass power 
spectrum compu ted from a n on-linear estimation of dark matter 
clustering dSmith et all 120031) . The two terms, <?k(x)> are me ge- 
ometric lens-efficiency, which depend on the redshift distribution 
of the sources, rik(x')' 



9k (x) 



X kU) /x(x') ' 



(13) 



Given a cosmological model, matter power spectrum, and red- 
shift distribution of the sources we can model the shear correlation 
functions. Bayesian model fitting techniques are then used to obtain 
the posterior probability on the model vector given the observed 
shear correlation functions. We discuss this further in Section[3~3l 



3.2 The tomographic weak lensing signal 

Based on the results presented in Section|2] the redshift distribution 
in each bin is taken to be the sum of the PDFs determined from 
the photometric redshift analysis of Hildebrandt et al.l d2012l) . We 
refer to the maximum posterior photometric redshift estimate as the 
'photometric redshift'. The histogram of photometric redshifts and 
the sum of the PDFs for each redshift bin are presented in Figure|5] 
Note that the summed PDFs extend to lower and higher redshifts 
then the photometric redshifts do, broadening the range below z = 
0.5 and above z — 1.3. The summed PDFs for the two redshift bins 
also overlap considerably with one another. The average redshift 
from the summed PDFs is 0.7 for the low redshift bin and 1.05 for 
the high redshift bin. For the photometric redshifts we find 0.69 and 

1.03 for the low and high redshift bins respectively. The average 
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Figure 5. Redshift distributions used in the weak lensing analysis. Low and 
high redshift bins correspond to z p = (0.5, 0.85] and z p = (0.85, 1.3] 
respectively. Smooth curves show the result of summing the photometric 
redshift probability distribution functions (PDFs) of all galaxies within the 
respective redshift bin, and the solid and dashed curves are used in the to- 
mographic analysis. The sum of the PDFs over the entire redshift range is 
given by the dot-dashed line which is used in the 2D lensing analysis. For 
comparison, the histograms show the redshift distribution obtained from the 
photometric redshifts. 



redshift for both bins taken together is found to be 0.87 from the 
summed PDFs and 0.84 from the photometric redshifts. 

We use ATHENA with an opening angle of 0.02 to measure 
the shear-shear correlation function. We have tested that the differ- 
ence to the shear-shear correlation function when using an opening 
angle of 0.02 compared to a brute force calculation is negligible, 
approximately 8 per cent of the size of the errors. The signal is first 
measured on each of the four wide mosaics: W1,W2,W3, and W4, 
applying the shear calibration described in Section [37X1 The cor- 
relation functions are then combined by calculating the weighted 
average. The weight for a given angular bin and wide mosaic is the 
sum of the inverse variance weight terms for each pair of galax- 
ies. We present and £_ for each redshift bin combination in 
Figure|6] The error bars correspond to the diagonal elements from 
the covariance matrix, discussed in more detail in Section |3,3, 11 
The lines are the theoretical prediction for a fiducial cosmologica l 
model using WMAP7 best-fitting results dKomatsu et all l201ll) . 
hence the following parameter vector is used: (f2 m = 0.271, 
a 8 = 0.78, h = 0.704, D. h = 0.0455, n s = 0.967, fi A = 0.729, 
Wo = —1). Descriptions of each parameter can be found in TableQ] 
To compute the t heoretical models we employ the halo-model of 
ISmith et all (2003) to estimate the non-li near matter powe r spec- 
trum and the analytical approximation of lEisenstein & Hul dl998h 
to estimate the transfer function. 

Emphasizing that no cosm ology-dependent syste matic tests 
were used to vet the catalogues (Hevmans et al, 2012b), Figure [6] 
demonstrates the robustness of the CFHTLenS catalogues. The to- 
mographic shear signal shows no evidence of a re dshift-dependent 
bias a s was seen in earlier CFHTLS data analyses (Kilbinger et al., 
2009). We discuss further tests of the redshift scaling of the shear 
in Section l3.3.3l 



3.3 Cosmology 

From the signal measured in Section 13.21 cosmological parame- 
ters are estimated using COSMOPMC. COSMOPMC is a freely 
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Table 1. Details of the model dependent cosmological parameters for each of the considered cosmologies. Parameter ranges denote hard priors. A flat 
distribution is used throughout the range. The bottom three parameters are constrained by WMAP7, and are required in order to deduce ag. 



Parameter flat ACDM curved ACDM description 



£2 m 


[0,1.0] 


[0, 1.2] 


Energy density of matter (baryons + dark matter). 




[0.2, 1.5] 


[0.2,1.5] 


Normalisation of the matter power spectrum. 


h 


[0.4,1.2] 


[0.4,1.2] 


The dimensionless Hubble parameter h = — — H 9 - , r . 


Q b 


[0,0.1] 


[0,0.1] 


Energy density of baryons. 


n s 


[0.7, 1.3] 


[0.7,1.3] 


Slope of the primordial matter power spectrum. 


U A 


i - n m 


[0,2] 


Energy density of dark energy. 


Wo 


-i 


-1 


Constant term in the dark energy equation of state, w (a) = wo 



T 


[0.04,0.20] 


[0.04,0.20] 


Reionisation optical depth. 


A 2 


[1.8,3.5] 


[1.8,3.5] 


Amplitude of curvature perturbations, units of 10 — 9 times the amplitude of density fluctuations. 


A sz 


[0.0,2.0] 


[0.0,2.0] 


Sunyaev-Zel'dovich template amplitude. 



availabl^EI Population Monte Carlo (PMC) code, which uses 
adaptive importance s ampling to explore the posterior likelihood 
(Kilbi nger et al , 2011 1. COSMOPMC documentation can be found 
in lKilbinger et alJ d201 lh : discussion of Bayesian evidence and ex- 
amples of its application to va rious cosmological data sets can be 
found in lKilbinger et alJ feOld) andlWraith etal] j2009h . The PMC 
method is detailed in ICappe et alJ (120070 . The non-l inear matter 
power spectrum is estimated using the halo-model of ISmith et alJ 
d2003h . The tr ansfer function is estima ted using the analytical ap- 
proximation of lEisenstein & Hu (1998). 

We explore two cosmologies: a flat ACDM universe and 
a curved ACDM universe. The model-dependent data vector 
used with COSMOPMC contains the following seven parameters: 
(fi m , ag, h, Sib, n s , SIa, wo)- Physical descriptions and priors are 
presented in Table [T] For the flat ACDM model we have a five- 
parameter fit, where we fix Q,a = 1 — fi m and too = — 1. For the 
curved ACDM model we have six free parameters as Qa is allowed 
to vary, while wo remains fixed. 

Three other cosmological data sets are used to provide com- 
plementary constraining power. Constraints from the cosmic mi- 
crowave ba ckground (CMB) are taken from the seven-year results 
of WMAP ( Kom atsu et"aill201ll . hereafter referred to as WMAP7). 
Baryon aco ustic oscillation (BAD) data is taken from the BOSS 
experiment dAnderson et al hereafter referred to as BOSS). 

The Hubble cons tant is constraine d with the results from the HST 
distance ladder (Ries s~et"all 1201 lL hereafter referred to as Rll). 
Following Rl 1, we use a Gaussian prior of mean value Ho — 0.738 
and sta ndard deviation a = .024. For more details of these data 
sets see Kilbinge r et alJ J2012t) . With WMAP7 the parameter set is 
expanded to include r, Asz, and A^, from which we deduce ag. 
Prior range s and brief description s are given in Table [T] For further 
details see Kom atsu et alJ feOllh and references therein. Through- 
out this section when stating parameter values we quote the 68.3 
per cent confidence level as the associated uncertainty with all other 
parameters marginalised over. 



Table 2. Constraints orthogonal to the Q m — erg degeneracy for a flat 
ACDM cosmology. Results are shown with and without highly non-linear 
scales which are potentially biased due to non-linear modelling and the ef- 
fects of baryons (see Sectionf4). 'All scales' refers to scales the correlation 
functions are measured on: 1 < 9 < 40 arcmin. We remove scales cor- 
responding to < 10 arcmin in the case labelled 'removed: < 10 
arcmin' . 



Data 






a 


tomography: 








all scales 




0.771 ±0.040 


0.553 ± 0.016 


removed: £_ 


< 10 arcmin 


0.776 ± 0.041 


0.556 ± 0.018 


2D Lensing: 
all scales 
removed: £_ 


< 10 arcmin 


0.785 ±0.036 
0.780 ± 0.043 


0.556 ± 0.018 
0.611 ± 0.015 



surveys constructed from the three-dimensional N - body n umer- 
ical lensing simulations of lHarnois- Dera ps et al.l d2012l) . The 
1024 3 particle simulations have a box size of 147.0 h _1 Mpc or 
231.1 h _1 Mpc, depending on the redshift of the simulation, and 
assume a flat A CDM cosmology para metrized by the best-fitting 
constraints from Ko matsu et alJ d2009l) . There are a total of 184 
fully independent lines of sight spanning 12.84 square degrees with 
a resolution of 0.2 arcmin sampled at 26 redshift slices between 
< z < 3. The two-point shear statistics match the theoretical 
predictions of the inp ut cosmology from 0.5 < 8 < 40 arcmin 
scales at all redshifts (Harnois-De raps et al .ll2012h . this sets the up- 
per ang ular limit for our tomographic analysis. See iHevmans et all 
d2012al) for a detailed discussion of covaria nce matrix estimation 
from the N-body simulations presented in Harnois-Deraps et"ai] 
d2012l) . including the required [Andersor] d2003l) correction that we 
apply to de-bias our estimate of the inverse covariance matrix used 
in the likelihood analysis that follows. 



3.3.1 Covariance Matrix 

In order to estimate a covariance matrix for our measured shear 
correlation functions in equation [8] we analyse mock CFHTLenS 

5 http://cosmopmc.info 



3.3.2 Flat ACDM 

We present marginalised two-dimensional likelihood constraints in 
the Q m — ag plane in Figure|7] The best constraint from weak lens- 
ing alone is for a combination of f2 m and ag, which parametrizes 
the degeneracy. We find ag (§f?) a = 0.771 ± 0.040 with a = 
0.553 ±0.016. 
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Figure 6. The filled circles with solid lines and rilled squares with dashed 
lines show the measured signal for £4. and §— respectively. Each panel 
shows the shear correlation functions for a unique pairing of redshift bins. 
The top, middle and bottom panels correspond to low redshift correlated 
with low redshift (low-low), low with high redshift (low-high), and high 
with high redshift (high-high). Error bars are the square-root of the diag- 
onal of the covariance matrix measured from mock catalo gues (see Sec- 
tion l3.3.H . Theoretical predictions for a fiducial (WMAP7, Koma tsu et all 
|201 Vt) cosmology are presented as lines; these are not the best-fitting mod- 
els. There are two negative data points for §_ in the top panel, their values 
are -2.3 X and -4.9 X 10~ 6 for scales 1.34 and 2.18 arcminutes 

respectively. 



When combining CFHTLenS with WMAP7, BOSS, and Rl 1 
data sets, we findf2 m = 0.2762±0.0074 and cr 8 = 0.802±0.013. 
The precision is ~20 times better than for CFHTLenS alone where 
we find n m = 0.27 ± 0.17 and a 8 = 0.67 ± 0.23. Constraints 
on the full set of parameters are presented in Table [3] We show 
the results for CFHTLenS tomography, CFHTLenS combined with 
WMAP7, BOSS and Rll, and, to assess the contribution of our 
data set to these constraints, we include results for WMAP7 com- 
bined with Rll and BOSS. The most valuable contribution from 
CFHTLenS is for f2 m , as, and Qb, where we improve the preci- 
sion of the constraints by an average factor of 1.5. 
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Figure 7. Marginalised parameter constraints (68.3, 95.5, and 99.7 per cent 
confidence levels) in the C m — ag plane for a flat ACDM model. Results 
are shown for CFHTLenS (blue), WMAP7 (green), CFHTLenS combined 
with WMAP7 (black), and CFHTLenS combined with WMAP7, BOSS and 
Rll (pink). 



For comparison we perform the analysis with a single redshift 
bin spanning the range of our 2-bin analysis, 0.5 < z p ^ 1.3. 
We refer to this as the 2D lensing case, in contrast to the tomo- 
graphic case where we split the galaxies into two redshift bins. Fig- 
ure[8]shows the marginalised parameter constraints in the fi m — as 
plane for both 2D lensing and tomography, and the two cases result 
in very similar constraints. For 2D lensing we find as {w^f)" — 
0.785 ± 0.036 and a = 0.556 ± 0.018, which is in agreement 
with what we find for tomography (Table [2}. When combining the 
2D lesning results from CFHTLenS with WMAP7, BOSS, and Rl 1 
data sets, we find fi m = 0.2774±0.0074 and a 8 = 0.810±0.013, 
which are nearly identical to those found for tomography (listed 
above and in Table [3). This level of agreement is also found for 
all other parameters when combining CFHTLenS with the other 
data sets. We note that all the parameter estimates agree within the 
68.3 per cent errors and the size of the error bars from 2D lensing 
are very similar to those found with tomography when combining 
CFHTLenS with WMAP7, BOSS, and Rl 1 data sets. 

For CFHTLenS alone the parameter estimates for 2D lensing 
and tomographic lensing agree with each other within their 68.3 
per cent uncertainties, however, the parameter constraints do not 
improve. With two broad overlapping redshift bins of average red- 
shift 0.7 and 1.05, there appears to be insufficient additional in- 
formation to tighten parameter constraints. Previous estimates of 
the improvement in constraints from weak lensing tomography use 
non-overlapping redshift bins, Gaussian covariance, and es timate 
errors using a Fisher matrix analysis (see ISimon et all 120041) . For 
two redshift bins with z < 3 and divided at z — 0.75 they find 
the error on individual parameters from tomography to be 88 per 
cent those from a 2D analysis. With our overlapping redshift bins 
and non-Gaussian covariance it is not surprising that this marginal 
improvement is significantly degraded. Additionally, our 2D lens- 
ing result is for redshifts 0.5 < z p < 1.3, this removes low red- 
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Table 3. Parameter constraints with 68.3 per cent confidence limits. The 
following parameters are deduced for CFHTLenS: and qo. When com- 
bining data sets the deduced parameters are: erg, Q^, and qo- The label 
CFHTLenS+Others refers to the combination of CFHTLenS, WMAP7, 
BOSS, and Rll. 



Parameter 


flat ACDM 


curved ACDM 


Data 




0.27 ±0.17 
0.288 ±0.010 
0.2762 ± 0.0074 


0.28 ±0.17 
0.285 ±0.014 
0.2736 ± 0.0085 


CFHTLenS 
WMAP7+BOSS+R1 1 
CFHTLenS+Others 


CT8 


0.67 ± 0.23 
0.828 ± 0.023 
0.802 ± 0.013 


0.69 ±0.29 
0.819 ±0.036 
0.795 ±0.013 


CFHTLenS 
WMAP7+BOSS+R1 1 
CFHTLenS+Others 




1 - On 
1 -fim 
l-Qm 


0.38 ±0.36 
0.717 ±0.019 
0.7312 ± 0.0094 


CFHTLenS 
WMAP7+BOSS+R1 1 
CFHTLenS+Others 


n K 







0.19 ±0.43 CFHTLenS 
-0.0020 ± 0.0061WMAP7+BOSS+R11 
-0.0042 ± 0.0040CFHTLenS+Others 


h 


0.84 ± 0.25 
0.692 ± 0.0088 
0.6971 ± 0.0081 


0.81 ±0.24 
0.694 ±0.012 
0.693 ±0.011 


CFHTLenS 
WMAP7+BOSS+R1 1 
CFHTLenS+Others 


n b 


0.030 ±0.029 0.031 ±0.030 
0.0471 ± 0.0012 0.0472 ± 0.0016 
0.04595 ± 0.00086 0.0470 ± 0.0015 


CFHTLenS 
WMAP7+BOSS+R1 1 
CFHTLenS+Others 


'/() 


-0.57 ±0.27 
-0.568 ±0.016 
-0.585 ±0.011 


-0.29 ±0.40 
-0.574 ±0.025 
-0.594 ±0.014 


CFHTLenS 
WMAP7+BOSS+R1 1 
CFHTLenS+Others 


n s 


0.93 ±0.17 
0.965 ±0.012 
0.960 ±0.011 


0.91 ±0.17 
0.969 ±0.014 
0.972 ±0.012 


CFHTLenS 
WMAP7+BOSS+R1 1 
CFHTLenS+Others 


T 


0.086 ±0.014 
0.081 ± 0.013 


0.086 ± 0.015 
0.085 ±0.015 


WMAP7+BOSS+R1 1 
CFHTLenS+Others 


A 2 


2.465 ±0.086 
2.429 ±0.081 


2.45 ±0.13 
2.361 ±0.094 


WMAP7+BOSS+R1 1 
CFHTLenS+Others 




0.97 ±0.62 
1.33 ±0.60 


1.35 ±0.61 
1.39 ±0.57 


WMAP7+BOSS+R1 1 
CFHTLenS+Others 



shift galaxies with small signal-to-noise improving constraints and 
weakening the gains from dividing the redshift range. We therefore 
expect a modest improvement at best. We find that our tomographic 
errors are on average 1 16 per cent of our 2D lensing errors. The fact 
that we find larger errors for tomography is surprising and warrants 
further discussion. 

To test our covariance matrices we perform the analysis again, 
replacing the measured shear correlation function with that pre- 
dicted from our model using a WMAP7 cosmology. In this case 
we find that the tomographic errors are 98.2 per cent of the 2D er- 
rors. Therefore the increase in the tomographic errors compared to 
the 2D errors is not an inherent product of the covariance matrices 
used. This result confirms that for the case of overlapping redshift 
bins and non-Gaussian covariance the expected improvement from 
tomography is marginal, at best. 

As discussed in detail in section [4] we have also analysed the 
data after removing small scales which could be affected by errors 
in the non-linear modelling of the matter power spectrum and bary- 
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Figure 8. Marginalised parameter constraints (68.3 per cent confidence 
level) in the C m — as plane for a flat ACDM cosmology. We compare 
the results for 2D lensing (blue) and 2-bin tomography (green). We com- 
bine CFHTLenS with WMAP7, BOSS, and Rl 1 . Results are shown for 2D 
lensing (black) and 2-bin tomography (pink). 



onic effects. When removing these scales (£_ < 10 arcmin) we see 
an improvement in the errors for tomography which are measured 
to be 104 per cent those of 2D lensing with the same scales re- 
mo ved. The remaining discrepancy could be due to several factors. 
The lSmith et~ai1 J2003h non-linear prescription could easily be bi- 
ased at the few per cent level. Residual errors in the redshift of 
galaxies or other per cent level systematics could be present. In ad- 
dition we expect some degradation of the tomography errors due 
to the bias correction of the inverse covariance matrix dAndersonl 
2003). The covariance is estimated from a finite number of mock 
catalogues (see Section [3. 3. U , since the tomographic covariance 
contains three times the number of elements as the 2D covariance, 
measuring it fro m the same number o f mock catalogues results in a 
noisier measure. lHartlapetai1 l l2007l) predict an erroneous increase 
in likelihood area of three per cent given our number of mock cat- 
alogues 184 and the size of the data vector for 2D lensing 16 and 
tomography 48. 

Finally we note that our two redshift bins are chosen based 
on concerns of intrinsic alignment contamination, as such, they are 
not optimised for constraining cosmology. With more carefully se- 
lected redshift bins it may be possible to overcome the issues dis- 
cussed above and obtain improved cosmological constraints. 



3.3.3 Redshift scaling of the cosmic shear signal 

Previous CFHTLS data were found to underestimate the shear sig- 
nal at high redshift necessitating additional calibr ation parameters 
when performing cosmological fits to the data (Kilbinger et al., 
2009). We demonstrate here that the CFHTLenS data have a red- 
shift dependent shear-signal which agrees with expectations from 
the modelled ACDM cosmology. 

The excellent agreement between the 2D and tomographic 
lensing results (Figure [8} suggests that the shear signal across our 
two redshift bins is scaling as expected. This is also observed in the 
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Figure 9. Marginalised parameter constraints (68.3 per cent confidence 
level) in the f2 m — erg plane for a flat ACDM cosmology. The results are 
shown for each combination of the two redshift bins. The low and high 
redshift bins correspond to 0.5 < z p ^ 0.85 and 0.85 < z p ^ 1.3 re- 
spectively. The excellent agreement shows that redshift scaling of the signal 
is consistent with the modelled ACDM cosmology. 



excellent agreement between the measured shear and the shear pre- 
diction based on a fiducial WMAP7 cosmology shown in Figure[6] 

The shear correlation function for each pair of tomographic 
redshift bins is analysed separately, corresponding to the shear cor- 
relation functions shown in each panel of Figure [6] In Figure [9] 
we present marginalised parameter constraints (68.3 per cent confi- 
dence level) in the fi m — og plane for each redshift bin combination. 
Since each contour is obtained from a sub-sample of the full data- 
set the degeneracy between the parameters is more pronounced and 
the area of the contours is larger than when analysing the full data 
set (Figure [7}. The agreement between the contours in Figure [9] is 
a convincing demonstration that the redshift scaling of the shear in 
the CFHTLenS data is consistent with expectations from the mod- 
elled ACDM cosmology. 

The power-law fits to the degenerate parameter constraints in 
Figure[9]for each case are erg (£§7)° = 0.820 ± 0.067, 0.753 ± 
0.053, and 0.753 ± 0.050 with a = 0.662 ± 0.020, 0.621 ± 
0.016, and 0.535±0.013 for the low-low, low-high, and high-high 
redshift bin pairings respectively. Note the expected evolution of a 
with redshift. 

We reiterate that the cosmological model-dependent verifi- 
cation of redshift scaling presented here is completely indepen- 
dent of the calibration of the data, and the rejection of bad fields, 
that were done with te sts which are not sensitive to cosmology 
l lHevmans et~aill2012rj) . 



3.3.4 Curved ACDM 

A curved ACDM cosmology is modelled, for the full details of 
parameters and priors used see Table [TJ We present constraints in 
the Q m — as and fi m — Qa plane in Figure [10] One-dimensional 
marginalised results when combining CFHTLenS with WMAP7, 
BOSS, and Rll are fi m = 0.2736 ± 0.0085, fi A = 0.7312 ± 



0.0094, fi K = -0.0042 ± 0.0040, and erg = 0.795 ± 0.013. The 
constraints on fi m and erg do not change significantly from the flat 
ACDM case. Parameter constraints for both models are presented 
in Tabled The addition of CFHTLenS to WMAP7, BOSS, and Rl 1 
is most helpful at constraining fi m , erg, fix, and Qa- The precision 
for these parameters improves, on average, by a factor of two. 

We again find excellent agreement with the 2D lensing analy- 
sis. When combining the 2D lensing of CFHTLenS with WMAP7, 
BOSS, and Rll data sets, we find fi m = 0.2766 ± 0.0082, 
fi A = 0.7273 ± 0.0089, fi K = -0.0035 ± 0.0035 and erg = 
0.804 ± 0.016. We do not show the complete details of our 2D 
lensing parameter estimations. However, we note that in all cases, 
either with CFHTLenS alone or combined with the other cosmo- 
logical probes, the 2D results agree with the tomographic results 
within the 68.3 per cent errors and the size of the error bars are 
similar for both cases. We again find that for CFHTLenS alone the 
individual parameter uncertainties from the tomographic analysis 
are 1 16 per cent those of the 2D lensing analysis. 



3.3.5 Constraining the deceleration parameter 

The deceleration parameter go parametrizes the change in the ex- 
pansion rate of the Universe. We calculate this as a deduced param- 
eter for both the flat and the curved ACDM models. The decelera- 
tion parameter depends on the energy density parameters 

a(t )a(t ) _ fi n 



a 2 (t ) 



fi A (curved ACDM), and 



3fi n 



1 (flat ACDM), 



(14) 



where the scale factor at present time is a(to) and derivatives with 
respect to time are denoted with a dot. For the flat case qo is simply 
a transformation of our results for the matter density parameter fi m . 
We present marginalised constraints for go in Figure QTJ The pink 
line is for the curved case where we find go = —0.29 ± 0.40, 
and the blue line is for the flat case where we find go = —0.57 ± 
0.27. Negative values indicate acceleration of the expansion of the 
Universe. Summing the posterior for go < tells us the confidence 
level at which we have measured an accelerating Universe. For the 
curved and flat models we find that go < at the 82 and 89 per 
cent confidence level, respect ively. 

ISchrabback et alj d2010h constrain go with a six-bin tomo- 
graphic analysis of the COSMOS-30 data. Besides having more 
tomographic bins the redshift range probed is also greater, extend- 
ing to z = 4. For a curved ACDM cosmology, holding fib and 
n s fixed and using a Gaussian prior on the Hubble constant of 
h = 0.72 ± 0.025, they find go < at 96 per cent confidence. 
If we do a similar analysis with fib and n 3 held fixed and us- 
ing the a Gaussian prior of h = 0.738 ± 0.024 (Rll), we find 
go < at 84 per cent confidence. The difference in constraints 
on the deceleration parameter can be understood as a result of 
the much larger values of the dark-energy density preferred by 



COSMOS-30 fi A = 0.97Z 



which lead to smaller values of go . 



Whereas the dark-energy density found here from CFHTLenS is 
fiA = 0.38 ± 0.36, resulting in larger values of go for CFHTLenS. 

With CFHTLenS alone we are not able to put a strong con- 
straint on the acceleration of the expansion of the Universe. With 
the addition of the other cosmological probes the entire pos- 
terior distribution of go is less than zero. For a curved model 
with CFHTLenS combined with the other probes, we find go = 
—0.594 ± 0.014 (see Table|3j. An accelerating Universe is unam- 
biguously detected. 
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Figure 10. Marginalised parameter constraints (68.3, 95.5, and 99.7 per 
cent confidence levels) for a curved ACDM cosmology. Results are 
shown forCFHTLenS (blue), WMAP7 (green), CFHTLenS combined with 
WMAP7 (black), and CFHTLenS combined with WMAP7, BOSS, and Rl 1 
(pink). Top panel: Constraints in the Q m — erg parameter space. Bottom 
panel: Constraints in the f2 m — S1a parameter space. 

4 IMPACT OF NON-LINEAR EFFECTS AND BARYONS 
ON THE TOMOGRAPHIC COSMOLOGICAL 
CONSTRAINTS 

We have presented cosmological parameter constraints from an 
analysis of the tomographic two-point shear correlation function 
£± W (equation ! Ill , incorporating the n on-linear dark matter only 
power spectrum from lSmith et al.l d2003l) as our theoretical model 
of Ps(k, z) in equation J 1 2b - Note that k — — 4^ . This halo-model 
prescription for the non-linear correction has been calibrated on nu- 
merical simulations and shown t o have an a ccuracy of 5 — 10 per 
cent over a wide range of scales (Eiflerll201 ]]). The N-body simula- 
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Figure 11. Marginalised constraints on the deceleration parameter using the 
CFHTLenS 2-bin tomographic weak lensing results. An accelerating uni- 
verse (go < 0) is found at the 82 per cent confidence level for a curved 
ACDM model (pink), and at the 89 per cent confidence level for a flat 
ACDM model (blue). 



tions used to estimate the covariance matrices used in this analysis 
suggest t hat the accuracy is even bette r than this for a WMAP5 cos- 
mology dHarnois-Deraps et al 1 12011) over the redshift range cov- 
ered in this analysis. While these comparisons give us confidence 
in our results, and suggest that any error from the non-linear correc- 
tion will be small in comparison to our statistical error, it is prudent 
to assess how errors in the non-linear correction will impact our 
results. 

A fully 3D weak lensi n g ana lysis of the CFHTLenS data is 
presented in iKitching et all «2012h . This power spectrum analy- 
sis allows for exact redshift dependent cuts in the wave-vector k, 
which can be motivated by either the comparison of the power 
spectrum measured from N-body simulations to the non-linear pre- 
scription, or the selection of linear scales where the non-linear cor- 
rection is neglig ible. For real-sp ace statistics, as used in this pa- 
per and in lKilbinger et alj d2012l) . it is not possible to make an un- 
ambiguous separation of scales. The two-point correlation function 
£± (0) is related to the underlying matter power spectrum Ps(k, z) 
through integrals over k and z, modulated by the lensing efficiency 
g(z) (equationll3t and Bessel functions Jo/4((-d) for £+/_ (6) (see 
equations 1111 and [121 . The measured tomographic two-point shear 
correlation function at a fixed scale 6 is therefore probing a 
range of k in the underlying matter power spectrum. In addition, 
owing to the different Bessel functions, is preferentially prob- 
ing much smaller k, and he nce larger physical scales, than £_ . 

iKilbinger et aljfc012h present an analysis of the 2D shear cor- 
relation function out to large angular scales < 350 arcmin. The 
consistent constraints obtained from the large quasi-linear regime 
8 > 53 a rcmin in comparis on to the full angular range, analysed 
using the lSmithetaT] d2003l) non-linear power spectrum, give us 
confidence that the accuracy of this correction is sufficient, falling 
within our statistical errors. 

Comparing the theoretical expectation of f^' (0) (equationll It 
for a WMAP7 cosmology, calculated using a non-linear and a lin- 
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Figure 12. Marginalised parameter constraints (68.3, 95.5, and 99.7 per 
cent confidence levels) in the f2 m — erg plane for a flat ACDM cosmology. 
The pink contours show the result when all eight scales are included, this is 
the same as the result for CFHTLenS shown in Figure|7] The blue contours 
show the result of removing highly non-linear scales, which are possibly 
biased due to the non-linear correction to the matter power spectrum or the 
effect of baryons. We remove the 5 smallest scales of corresponding to 
< 10 arcmin. The contours are only slightly different, indicating that we 
are not sensitive to these effects given the level of precision of our results. 

ear power spectrum, we determine the angular scale below which 
the non-linear and linear models differ in amplitude by greater than 
10 per cent. For £+, this quasi-linear limit ranges from 10 — 14 
arcmin for the three different tomographic combinations (the low- 
est redshift bin requiring the largest 6 cut). For £_ , the quasi-linear 
limit ranges from 100 — 140 arcmin. In this analysis, we limit our 
angular range to scales with 6 < 40 arcmin where we can ac- 
curately assess a covariance matrix from the lensing simulations 
dHarnois-Deraps et all l2012h . We are therefore unable to follow 
Kilbingeret al. (2012) by limiting our real-space analysis to this 
quasi-linear regime as we do not probe sufficiently large angular 
scales. We can however make an assessment of how an error on the 
non-linear correction would impact our results. We first compare 
the WMAP7 theoretical expectation of C±'(^) calculated using a 
non-linear correction boosted by 7 per cent, with a model calcu- 
lated with the non-linear correction decreased by 7 per cent. Note 
that we chose the valu e of 7 per ce nt from the average error over 
the range of k tested in lEiflerl J20T ]]). We find that these two limits 
on the non-linear correction cause at least a 10 per cent change in 
the amplitude of £±'(0) for scales 9 < 1 arcmin for and 9 < 10 
arcmin for Applying these cuts in angular scale corresponds 
to removing the first 5 angular scales for £_ for each tomographic 
bin shown in Figure [6] All scales remain in the analysis. With 
these scales removed any remaining uncertainty due to the non- 
linear modelling is well within our statistical error. 

Figure[T2]compares cosmological parameter constraints in the 
S^m — o"8 plane for this limited number of scales in comparison 
to the full data set analysed in Section [3] The removal of small 
scales results in a slight change to the degeneracy of the parameters. 
This test gives us confidence that the non-linear correction used is 



sufficiently accurate given the statistical error of the survey. This is 
unlikely to be true for future surveys, where the increased statistical 
accuracy will require b etter knowl edge of the non-linear correction 
to the power spectrum (Eifle tolll) . 

Finally we turn to the impact of baryons on our results. In our 
analysis we assume the underlying matter power spectrum is suffi- 
ciently well represented by the non-linear dark matter only power 
spectrum, neglecting the role of baryons. The impact of baryons 
on the power spectrum is sensitive to the baryonic feedback model 
used. The refore, the magni t ude of the impact of baryons remains 
uncertain. Semboloni et al. I fcOllh present an analysis of cosmo- 
logical hydrodynamic simulations to quantify the effect of baryon 
physics on the weak gravitational lensing shear signal using a range 
of different baryonic feedback models. Their work suggests that 
a conservative weak lensing analysis should be limited to those 
scales where k < 1.5hMpc _1 . We implement su ch a conserva- 
tive sc heme in the 3D power spectrum analysis of iKitching et alj 
d2012h . As discussed above, our real-space analysis mixes k and z 
scal es, leaving us unable to p erform a similarly clear test here. 

ISemboloni et al.l J20 1 lh present a comparison of f ± ' (0) mea- 
sured for both the cosmological hydrodynamic simulations and a 
dark matter only simulation for different redshifts, which we use 
to judge the level of error we should expect baryons to introduce. 
Assuming the realistic AGN feedback model, and considering the 
scales used in the conservative analysis of Figure [T2l (9 1.34 
arcmin and 9 ^ 15.4 arcmin for and respectively), we ex- 
pect baryons to cause a decrease to the modelled signal of less than 
ten per cent. The fact that we see very little difference between the 
conservative and full analysis presented in Figure [L2ldemonstrates 
that the underlying matter power spectrum is indeed sufficiently 
well represented by the non-linear dark matter only power spec- 
trum for our statistical accuracy. It also indicates that the impact of 
baryons on the non-linear dark matter only power spectrum is un - 
likely to be larger than that predicted by Sem boloni et al.l d201ll) . 
However, baryonic effects will have to be carefully considered for 
the next generation of weak len sing surveys that will ha ve signif- 
icantly smaller statistical errors. ISemboloni eTail d2012l) show the 
importance of baryonic effects on three-point shear statistics and 
propose a modification to the modelling of the non-linear matter 
power spectrum to account for these effects. 



5 CONCLUSION 

The most important result of this study is that the sum of the pho- 
tometric redshift probability distribution functions (PDF) within a 
redshift bin provides an accurate measure of the true redshift dis- 
tribution of those galaxies; accounting for the scatter due to catas- 
trophic as well as statistical errors. To demonstrate the accuracy 
of the PDFs we have compared the summed PDFs with the red- 
shift distribution predicted by spectroscopic redshifts, resampled 
COSMOS-30 redshifts, and predictions from a redshift contamina- 
tion analysis using the angular correlation function. We find excel- 
lent agreement for the redshift range z p < 1.3. T his result indicates 
that th e priors and spectral templates used in [rlildebran dt et alj 
d2012l) to derive the photometric redshifts provide an accurate and 
complete description of the galaxies at z p < 1.3. This also moti- 
vates our use of the summed PDF as a measure of the redshift distri- 
butions in our tomographic weak lensing analysis. Furthermore, the 
proven accuracy of the summed PDFs provides a reliable method 
for estimating the source redshift distribution in future weak lens- 
ing studies. 
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We have performed a cosmological analysis of the CFHTLenS 
data on angular scales 1 < 9 < 40 arcmin, using two broad 
redshift bins, 0.5 < z p < 0.85 and 0.85 < z p < 1.3, that 
are not significantly affected by the intrinsic alignment of galaxy 
shapes. We model two cosmologies; flat and curved ACDM. Due 
to complementary degeneracies our results add valuable constrain- 
ing power w hen combined with th ose from the cosmic microwave 
backgroundj Komatsu et all 201 1[ WMAP7), baryon acoustic os- 
cillations ( Ander son et all 20121 BOSS), and a prior on the Hubble 
parameter l lRiess et al. , 201 ll Rl 1). The addition of our weak lens- 
ing results to these other cosmological probes increases the preci- 
sion of individual marginalised parameter constraints by an average 
factor of 1.5 - 2. 

For a flat ACDM model the joint parameter constraints for 
CFHTLenS, WMAP7, BOSS, and Rl 1 are Q m = 0.2762±0.0074 
and a s = 0.802 ± 0.013. For a curved ACDM model, combin- 
ing the same data sets, we find fi m = 0.2736 ± 0.0085, Qa = 
0.7312 ± 0.0094, Q K = -0.0042 ± 0.0040, and a 8 = 0.795 ± 
0.013. Full details of our parameter estimates for both cosmologies 
are presented in Table [3] Our results are consistent with those pre- 
sented in other studies of the CFHTLenS data: a 2D lensing analy- 
sis probing much larger scales where linear theory provides a more 
accur ate model to the matter power spectrum (Kilbinger et al., 

HH3); 

and a fine-binned tomographic analysis with six red shift 
bins accounting for intrinsic alignments (IHevmans et al ., 2012a). 

We compare the tomographic constraints with those from a 2D 
lensing analysis spanning the same range of redshift 0.5 < z p ^ 
1.3. We find the two analyses to be completely consistent with all 
parameter estimates agreeing within their 68.3 per cent confidence 
levels. We note that the uncertainties on individual parameters from 
tomography are on average 116 per cent larger than the uncertain- 
ties from 2D lensing. This statistic is 98.2 per cent if we replace our 
data vectors with a fiducial model, indicating that our covariance 
matrices do show a slight improvement for tomography. We argue 
that our non-Guassian covariance and broad overlapping redshift 
bins degrade the modest improvement ( ~88 per cent) exp ected 
from idealised Fisher matrix calculations (Sim on et al 1 120041) . We 
identify small scales as being responsible for 12 per cent of the 
increase from 98.2 per cent. These scales could be biased due to 
uncertainties in the modelling of the non-linearities in the matter 
power spectrum and baryonic effects. Although small scales have 
inflated our uncertainties from tomography we show in Section [4] 
that they do not significantly affect our results. 

Previous analyses of CFHTLS data were hindered by a strong 
redshift dependent bias in the weak lensing shear, necessitating ad- 
ditional nuisance param eters when analysing the tomographic shear 
(Kilbi nger et"ail [2009). We demonstrate that the redshift scaling of 
the CFHTLenS cosmic shear signal agrees with expectations from 
the modelled ACDM cosmology. The strongest test of this is pre- 
sented in Figure [9] which shows the agreement of cosmological 
constraints measured for each combination of redshift bins in the 
(J m — as plane for a flat ACDM cosmology. This demonstrates the 
effectiveness of the c osmology-independent tests of residual sys- 
tematics presented in Hev mans et all J2012bl) . including the rejec- 
tion of 25 per cent of the MEGACAM pointings which failed to 
pass these tests. Note a lso that the shear c alibration performed on 
numerical simulations (Mil ler et al 1 12012b was completed before 
any cosmological analysis was performed on the data. The two-bin 
analysis presented here is sensitive to redshift dependent cosmol- 
ogy without introducing additional parameters to model intrinsic 
alignments, as such it is an excellent final test of the CFHTLenS 
data product. 
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